Single-photon nonlinearities in two-mode optomechanics 
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We present a detailed theoretical analysis of a weakly driven multimode optomechanical system, 
in which two optical modes are strongly and near-resonantly coupled to a single mechanical mode 
via a three-wave mixing interaction. We calculate one- and two-time intensity correlations of the 
two optical fields and compare them to analogous correlations in atom-cavity systems. Nonclassical 
photon correlations arise when the optomechanical coupling g exceeds the cavity decay rate k, and we 
discuss signatures of one- and two-photon resonances as well as quantum interference. We also find 
a long-lived correlation that decays slowly with the mechanical decay rate 7, reflecting the heralded 
preparation of a single phonon state after detection of a photon. Our results provide insight into the 
quantum regime of multimode optomechanics, with potential applications for quantum information 
processing with photons and phonons. 

PACS numbers: 42.50.Wk, 42.50.Lc, 07. 10. Cm 



I. INTRODUCTION 

Optomechanical systems (OMSs) involve the interac- 
tion between optical and mechanical modes arising from 
radiation pressure force, canonically in an optical cav- 
ity with a movable mirror [I^Sj. Recent progress in op- 
tomechanical (OM) cooling techniques has been rapid ,4r 
[TT| . and experiments have now demonstrated cooling to 
the mechanical ground state [T^HTi] . OM induced trans- 
parency |151 I16j . and coherent photon-phonon conver- 
sion [T7Ul9j . These developments have attracted signif- 
icant interest, and motivated proposals for applications 
exploiting OM interactions at the quantum level, ranging 
from quantum transducers [20,-23^ and mechanical stor- 
age of light [24H26] to single-photon sources [27] and OM 
quantum information processing [281 129] . Significant ad- 
vantages of OM platforms for these applications are the 
possibility for mass production and on-chip integration 
using nanofabrication technologies, wide tuneability and 
the versatility of mechanical oscillators to couple to a 
wide variety of quantum systems |30) . 

The force exerted by a single photon on a macroscopic 
object is typically weak; consequently, experiments have 
so far focused on the regime of strong optical driving, 
where the OM interaction is strongly enhanced but ef- 
fectively linear [3T1 [35] However, recent progress in the 
design of nanoscale OMSs [TT, '33'-'3F and OM experi- 
ments in cold atomic systems [36t |37j suggests that the 
regime of single-photon strong coupling, where the OM 
coupling strength g exceeds the optical cavity decay rate 
K, is within reach of the next generation of OM experi- 
ments. In this regime, the inherently nonlinear OM in- 
teraction is significant at the level of single photons and 
phonons JT, 38 40J. For example, the presence of a sin- 
gle photon can — via the mechanical mode — strongly in- 
fiuence or even block the transmission of a second photon. 



leading to photon blockade. This single-photon nonlin- 
earity was recently analyzed for canonical OMSs consist- 
ing of a single optical mode coupled to a mechanical mode 
[171 mi [55]. However, with a single optical mode, the OM 
coupling is highly off-resonant, leading to a suppression 
of effective photon-photon interactions by the large me- 
chanical frequency iOm ^ g |27j . 

In this paper we develop a quantum theory of a weakly 
driven two-mode OMS [281 l43ti45] in which two optical 
modes are coupled to a mechanical mode. The key advan- 
tage of this approach is that photons in the two optical 
modes can be resonantly exchanged by absorbing or emit- 
ting a phonon via three-mode mixing. We extend our ear- 
lier results [28 , where we discussed possible applications 
of resonant optomechanics such as single-photon sources 
and quantum gates, by exploring one-time and two-time 
photon correlations of both optical modes. Specifically, 
we find that the photon-photon correlation function of 
the undriven optical mode exhibits delayed bunching for 
long delay times, arising from a heralded single mechani- 
cal excitation after detection of a photon in the undriven 
mode. Finally, we compare the two-mode OMS to the 
canonical atomic cavity QED system with a similar low- 
energy level spectrum [46l [47] . Despite several similar- 
ities we find that, in stark contrast to the atom-cavity 
system, the OMS studied here does not exhibit nonclas- 
sical correlations unless the strict strong coupling condi- 
tion 5 > K is met. Our results serve as a guideline for OM 
experiments nearing the regime of single-photon nonlin- 
earity, and for potential quantum information processing 
applications with photons and phonons. 

The remainder of the paper is organized as follows. In 
Sec.[ll| we introduce the system and details of the model. 
In Sec. [III[ we calculate the equal-time intensity correla- 
tion functions of both transmitted and reflected photons, 
and discuss signatures of nonclassical photon statistics. 



In Sec. IV we investigate two-time correlation functions 
of the transmitted photons, and discuss delayed coinci- 
dence correlations that indicate the heralded preparation 
of a single phonon state. Finally, we provide a brief out- 
look on the feasibility of strong OM coupling in Sec. [Vj 
and conclude in Sec. |VI| with a summary of our results. 
The Appendix contains details of our analytic model used 
to derive several results discussed in the paper. 



II. MULTIMODE OPTOMECHANICS 

We consider the setup shown schematically in Fig.[TJa), 
consisting of two optical cavities separated by a semi- 
transparent mirror. The cavity modes are coupled by 
photons tunneling through the fixed mirror in the mid- 
dle, and the mode on the right couples to the motion of 
a vibrating endmirror through radiation pressure. The 
Hamiltonian describing the system is {h = 1) 



Ho = a;o(clci + c|c2) - J(c|c2 + c\c2) 

+ uj^b^b - g{b^ + b)c\c2 + Hdr{t), 



(1) 



where ci_2 are annihilation operators for the two cavity 
modes, which we assume to be degenerate with frequency 
wq, and J is the photon tunneling amplitude through the 
central mirror. The motion of the endmirror on the right 
is described by a single mechanical mode with annihila- 
tion operator b and frequency ujm, and the parametric 
coupling strength g corresponds to the shift of the cavity 
frequency due to a single mechanical phonon. Finally, 
Hdrit) = J2i=i2 (f^iCiC*'^^* -|- h.c.) describes two coher- 
ent driving fields of amplitudes fij and frequency wl, 
which are applied to the left and right cavities. 

We are interested in a three-mode resonant interac- 
tion in which the two optical modes exchange a pho- 
ton by absorbing or emitting a phonon in the mechani- 
cal mode. We begin by diagonalizing the optical part of 
Hq in the first line of Eq. ([l]) in terms of the symmetric 
and antisymmetric combinations of the optical modes, 
Cs = ^(ci -I- C2) and Ca = ^(ci - C2), with eigenfre- 
quencies LLia,s — ujq ± J . In the frame rotating at the 
laser frequency cj^ we obtain 



H' = - Acica - (A -t- 2J)clcs ^ 
+ f (ck + cK) (5 + 6t) 



LUrnbh 



(2) 



where A = ojl — LUa is the laser detuning from the Ca 
mode, and ils,a = {^1 i il.2)/V^- Next, we focus on the 
case of three-mode resonance, a;„i — 2 J, and assume that 
3> g, |A|, n^. This allows us to make a rotating wave 
approximation with respect to the remaining large fre- 
quency scale Wm, and in the frame defined by the trans- 
formation U = exp[—iuJm.t{b^b— cjcs)], the Hamiltonian 
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FIG. 1. (a) Optomechanical system consisting of two tunnel- 
coupled optical cavity modes ci and C2, and a mechanical 
oscillator b coupled to one of the cavity modes by radia- 
tion pressure. The coupled optical modes are diagonalized 
in terms of symmetric and antisymmetric modes, Cs and Ca- 
(b) Level diagram showing the relevant zero-, one- and two- 
photon states at zero temperature and under the three-mode 
resonance condition cJs — cOa = oJm- States are labeled by 
InaKsW;,) denoting the number Ua (ris) of antisymmetric (sym- 
metric) photons and the number of phonons ni,. The optome- 
chanical coupling g splits the degeneracy of states \nansni,) 
and \na — l,ns + l,nb + 1). 



H' simplifies to 



H = ~A{clca + clcs) + |(clc,& + b^clca) + l^a(cl -f- Ca) . 



(3) 



This is the starting point for our analysis discussed below. 
Note that the assumptions made in deriving Eq. (13) are 
fulfilled in most experimental systems of interest [48l [49] , 
and H represents a generic description for resonant two- 
mode optomechanics [131 [HI [SOI [M] • 

The nonlinear terms proportional to g in Eq. ([s]) de- 
scribe coherent photon exchange between the two op- 
tical modes mediated by absorption or emission of a 
phonon. The resulting low energy level diagram is shown 
in Fig. [l]jb), where \nansni,) represents a state with 
and Hs photons in the Ca and Cs modes, and phonons 
in the mechanical mode. In the absence of a drive we 
diagonalize H in this few-photon subspace, yielding the 
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eigenstates 
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(|100)±|011)), 
(|200) ±\/3|lll) 
\/2|200) - |022) 



721022)) 



(4) 
(5) 

(6) 
(7) 



Note that in the diagonal basis the weak driving field 
couples all states with photon number differing by one. 
In the following sections we use this low energy Hilbert 
space to understand photon correlations in the system. 

In addition to the coherent evolution modeled by the 
Hamiltonian H, we describe optical and mechanical dis- 
sipation using a master equation for the system density 
operator p, 



P - 



7 



K'D[Ca]p + k'D[Cs\p 



{Nt^ + l)V[h]p+lN,^V[b'^]p, 



(8) 
(9) 



2 

where H is given by Eq. (|3|, 2k and 7 are energy 
decay rates for the optical and mechanical modes, re- 
spectively, A'th is the thermal phonon population and 
'D[o\p = Ibpo^ — bp — po^ o. Below we study nonlinear 
effects at the level of single photons, both numerically 
and analytically, by solving Eq. ([s]) approximately in the 
limit of weak optical driving, f2 = ^ k. 



III. EQUAL-TIME CORRELATIONS 

A. Average transmission and reflection 

Before focussing on photon-photon correlations, we 
first study the average transmission through the cavity, 
which is proportional to the mean intracavity photon 
number. In Fig. [2]ja) and (b) we show the intracavity 
photon number of the two optical modes (dashed green 
curves) , 



{c]ci 



(10) 



where i = a, s, and angle brackets denote the steady 
state average. At A/g — ±i, both transmission curves 
exhibit a maximum, indicating that the driving field is in 
resonance with an eigenmode of the system. The position 
of these peaks can be understood from the level diagram 
shown in Fig. [ijb), which at finite g shows a splitting of 
the lowest photonic states into a doublet, |1±) = (|100)± 
|011))/y2.^ 

In addition to the transmission, we plot the mean re- 
flected photon number in Fig. [2f^c). As discussed below, 
the reflected photon statistics can also exhibit signatures 
of nonlinearity. We calculate properties of the reflected 
light using the cinni Illicit ion opeicitor cn — Ca ~\~ i-, ob- 
tained from standard input-output relations for a sym- 
metric two-sided cavity. The dashed curve in Fig. [2jc) 
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FIG. 2. (a) Normalized average photon number (dashed 
green line) and photon-photon correlation function (solid blue 
line) for driven mode Ca as a function of laser detuning at zero 
temperature. The average photon number is normalized by 
no = {Q./nf'. (b) Same as in (a) for the undriven mode Cs, 
and (c) for the reflected field cr. In all plots we took g/K = 20 
and "f / n = 0.2. Dots mark features seen at values of detuning 
A/3 = (A), ^ (B), I (C) and ^ (D). Bottom panels A-D 
illustrate the origin of these features as explained in the text. 
Supressions and enhancements of specific levels is indicated 
by a red X and a green O on top of them, respectively. 



shows the reflected photon number, nji 



At 



A/ (7 = the average reflection has a minimum where 
the average transmission has a maximum. Note that in 
contrast to a single cavity, even on resonance the trans- 
mission probability is less than unity and the reflection 
probability remains finite. 



B. Intensity correlations 

To characterize nonclassical photon statistics in the 
light transmitted though the OMS we study the equal- 



4 



time photon-photon correlation functions, 



where A denote steady state amphtudes. We obtain 



(2). . _ jclcUc,) 



(11) 



where aU operators are evaluated at the same time and 
i = a, s, R. A normalized correlation of g-j (0) < 1 i^" 

(2) 

dicates photon anti-bunching, and the limit gl^ (0) — 
corresponds to the complete photon blockade regime in 
which two photons never occupy the cavity at the same 
time. The solid curves in Fig. 2 show g'^aiO), gsa\o) 

(2) . — 

and 5/^^(0) as a function of the laser detuning and in 
the limit of weak driving fl/n <C 1. The most pro- 
nounced features of these correlation functions occur at 
\A\/g = 0, ^ and as marked by dots A, B, C and 
D, respectively. As we explain in detail in the follow- 
ing analysis, we find that the photon bunching at A and 
anti-bunching at B are the result of destructive quantum 
interference, while the features at points C and D arise 
from one- and two-photon resonances. 

To gain insight into the two photon correlation func- 
tions shown in Fig. [2] , we develop an approximate an- 
alytic model for the system by considering only the six 
levels shown in Fig. llTb). Assuming that the system is 
initially prepared in^OOO), these are the only levels sig- 
nificantly populated by weakly driving the Ca mode. We 
make the ansatz 1451 



= Aooo|000)+^ioo|100)+Aoii|011) 
+ A200I2OO) + Anillll) + A022IO22), 



(12) 



and describe the dynamics by evolving {ip) under the 
action of the non-Hermitian Hamiltonian, H = H — 
i [KcjjCa + Kclcg + 2 ^^^] ■ This approach allows us to 
evaluate intensities up to order ft^ and two-point correla- 
tion up to order 17*, since the neglected quantum jumps 
lead to higher order corrections. By neglecting the typi- 
cally small mechanical decay rate 7 ^ k, the amplitudes 
in Eq. ( 12 ) then satisfy 



^000 
^100 
ion 

^200 
Am 

A()22 



0, 



-i|^ioo ~ HAoii, 



(13) 
(14) 
(15) 

(16) 



-i — Aiii - iV2QAioo - 2kA2oo, 
v2 

-i^A20Q - igAo22 - ii^Aon - 2k^iii, (17) 



2kA 



022, 



(18) 



where k — k — lA. It is straightforward to solve Eqs. 
([l3|l8) for the steady state amplitudes (see Appendix). 
To lowest order in Q/k the mean occupation numbers 
are fia = l^iooP, = |ioiiP and ur = \Aioq + ifJ/Kp, 
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(19) 



(20) 



(21) 



where i?i<-(a;) = [i^^ -|- (A - a;)^] [K^ + {A + uj)^] and 
no ~ [Q./tif'. From the factors Rk{oj) in the denomi- 
nators (numerators) in these expressions, we obtain the 
positions of the resonances (antiresonances) in the aver- 
age intracavity photon numbers, in excellent agreement 
with the numerical results shown in Fig. [2] Our six-level 
model also provides the equal-time correlations (see Ap- 
pendix) , 



ffi^Ho) 



(±) R^ (f ) 



Vx/8 

R,{0)R, [^g 

R^ .9) ' 

Rk (f ) RwK^/g'-' 



2 ?T - 



(22) 



(23) 



(24) 



[R^/ADY 

Again, these expressions are in agreement with the fea- 
tures seen in the numerical results in Fig. [2j The po- 
sitions of maxima and minima are seen directly by the 
arguments of the factors (w). Note that we assumed 



qJk <C 1 to obtain the simplified expressions in Eqs. ( 22 
24 1, but we retained the shift of order g{K/g)'^ in the 
argument in Eq. ( 24 1 because this shift is larger than the 



width of the antiresonance. 

We now discuss each feature in Fig. |2] in terms of 
our six-level model together with the diagonal basis in 
Eqs. pfzl) . First, at detuning A = (point A in Fig. [2]) 

we see gda{0) > 1, indicating bunching. This is due 
to destructive interference that suppresses the popula- 
tion in 1 100) (panel A in Fig. [2]), and can be under- 
stood as the system being driven into a dark state, 
\d) (X g|000) — iljOll), similar to electromagnetically in- 
duced transparency (EIT) [TSl [S5]. In the dark state, 
|011) remains populated, allowing transitions to |111) 
which in turn is strongly coupled to |200). The net re- 
sult is a relative suppression of the probability to have 
one photon compared to two photons in the driven mode, 
leading to bunching at A = 0. Second, at detuning 
A = g/V& (point B), mode Ca shows anti-bunching due to 
a suppressed two-photon probability. Again, this is due 
to destructive interference, or the presence of a dark state 
in which |200) remains unpopulated (panel B). Third, 
at detuning A = | (point C), all modes show anti- 
bunching. This is due to a one-photon resonant transi- 
tion |0) |1+) (panel C). Finally, at detuning A = ^g 



5 



(point D), both Ca and Cs show bunching due to a two- 
photon resonant transition |0) |2+) (panel D). 



C. Absence of two-photon resonance at A = 

At first glance, the level diagram in Fig.[l]together with 
bunching in Fig. [2](a) suggest a two-photon resonance at 
zero detuning A = 0, where the energy of the eigen- 
state |2o) is equal to the energy of two drive photons. 
However, as discussed above, the bunching at A = 
arises entirely from the suppression of a one-photon pop- 
ulation; further, we find that the expected two-photon 
resonance is cancelled by interference. This can be seen 
from a second order perturbative calculation of the two- 
photon Rabi frequency JIq ^ for the transition |0) — >■ |2o). 
The two-photon state |2o) can be populated by the drive 
ifdr = ^(c^ + Ca) from state |0) via two intermediate one- 
photon eigenstates, |1±) given by Eq. ([s]), with energies 
LOi^ — — A± I in the rotating frame. The resulting Rabi 
frequency is 



(2) 
0,2o 



E 



(2o|i?dr|n)(n|ffdr|0) 



(25) 



n=l_ 



which vanishes at A = as a consequence of destruc- 
tive interference between the two amplitudes. The ex- 
act cancellation is lifted by including finite dissipation 
and the full spectrum; nonetheless this simple argument 
shows that the expected two-photon resonance at A = 
is strongly suppressed. 

Further evidence of the absence of a two-photon reso- 
nance at A = is the lack of bunching in the undriven 
mode in Fig.[2jb). If there were a two-photon resonance, 
one would expect that bunching should also occur in the 
undriven mode, since the state |2o) involves both Ca and 

Cs modes. This is indeed the case at detuning A = 
(see point D in Fig. [2| , where both modes show bunch- 
ing as a result of two-photon resonance. In contrast, we 
see no bunching in the undriven mode at A = 0. This 
supports our conclusion that the observed bunching at 
A = arises from suppression of population in |100) due 
to interference, as discussed in Section [lII B[ and not from 
two-photon resonance. As discussed above, this interfer- 
ence does not suppress population in |011), so we do not 
expect bunching in the Cg mode from this effect. 

Finally, to confirm our intuitive picture we plot the 
intensity correlation function, 5tot(0) — ("-tot ('^-tot — 
l))/(ritot)^5 of the total photon number, n^ot = 'T-a + ^sj 
in the coupled OM system in Fig. [3j The probability to 
find one photon in the combined cavity is maximal at 
A/g = ±i due to one- photon resonance. Similarly, we 
observe antibunching at point C and bunching at point 
D, due to interference and two-photon resonance respec- 



tively, as discussed in Section III B However, we find nei- 
ther bunching nor antibunching at A = 0, demonstrating 
the absence of a two-photon resonance despite the fact 
that |2o) lies at the twice the drive frequency. 




FIG. 3. Average number and intensity correlation of total 
photon number in the coupled OM system. Parameters are 
the same as in Fig. [2] and dots mark the same detunings. 
One- and two-photon resonances are seen at C and D, but 
we see no bunching in the total photon number at A = 
(point A). This reflects the lack of t wo- photon resonance due 
to destructive interference (see Eq. (25l). 



D. Finite temperature 



So far in our analysis we have focused on the case where 
the mechanical system is prepared in its vibrational 
ground state, |0m). This condition can be achieved using 
high frequency resonators operated at cryogenic temper- 
atures [H], and in the limit of weak driving Q,/k <^ 1 
such that optical heating of the mechanical mode can be 
neglected. In the following we extend our analytic treat- 
ment to the case of finite temperature, and show that 
many of the nonclassical features are robust even in the 
presence of small but finite thermal occupation of the 
mechanical mode. 

To generalize our previous results we now consider a 
finite equilibrium occupation number iVth > of the me- 
chanical mode, but still assume that 7(A'th -I- 1) ^ n,g- 
Within this approximation we proceed as above, and 
make a similar six-level ansatz as in Eq. (12 1 for each 
phonon number n. 



IV-n) = Ao,o,„|0,0,n) + Ai,o,„|l,0,n) 

+ Aoa,„+i|0,l,n + l) + A2,o,n|2,0,n) (26) 
+ Ai4,„+i|l, 1, n + 1) + Ao,2,n+2|0, 2, n + 2), 



where \'ipn) includes states up to two photons that are 
connected by the weak drive and coupling g, starting 
from the state |00n). As shown in Fig. 4] the coupling 
between the states within each six-level subspace depends 
explicitly on the phonon number n. Following the same 
approach as above, the amplitudes in Eq. (261 evolve 



6 



according to 

^0,0,n — 0, 



(a) 



gVa + 4n 



(27) 



.9 



(28) 

^0,1, «+i = + l^i^o.n - kAo,i,„+i, (29) 



/n + 1 



(30) 



2 -^V 2 

— iilAQin^l — 2kAiin+l, 



(31) 



A, 



0,2,n+2 — —T-g 



n + 2 



^ 2KAo_2,n+2- (32) 



We solve for the steady state amplitudes within each sub- 
space n and average the result over the initial thermal 
phonon distribution, assuming no coupling between sub- 
spaces due to the small phonon relaxation rate. We ob- 
tain the average photon numbers 

ria ^^Cn\Aifi^„f, fis ^^Cn\Ao,l,n+lf, (33) 



where C„ = £-''^'^""(1 - e"'^'^"-) and = keT. Sim- 
ilarly, the g^- (0) functions are given by 



siaHO) -2^C„I^AnlV^a, 

n 

gi^\0) =2j2Cn\Ao,2,n+2\Vf 



(34) 
(35) 



We provide the expressions for the steady state ampli- 
tudes ^2,0 n and 2,n+2 in the Appendix. 

In Fig. 4 we plot the correlation functions, g^aa (0) and 

(2) — 

gss (0) for different thermal phonon numbers, iVth- Solid 
lines were calculated from the above analytic approach 
with 7 — > 0, and we find excellent agreement with the 
full numerical results including small but finite 7 (dots, 
shown only for thermal occupation A^th — 2). We see 
that the zero temperature features such as antibunching 
survive at finite temperature for sufficiently strong cou- 
pling |28j . In the insets we plot the minimum antibunch- 
ing as a function of thermal occupation number for sev- 
eral ratios g/n. In addition, for detunings | A| > g/2, a se- 
ries of new resonances appear in the correlation functions, 
and for small but finite occupation numbers we find new 
antibunching features that are absent for A'th = 0. These 
new features can be understood from the n-dependent 
splitting of the one- and two-photon manifolds as indi- 
cated in Fig. |4]ja). For higher temperatures the individ- 
ual resonances start to overlap, and we observe an overall 
increase over a broad region of large positive and nega- 
tive detunings due to the cumulative effect of different 
phonon numbers. 
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FIG. 4. Correlation functions for finite temperature, 
(a) Level diagram showing the six states populated by the 
drive from level |0, 0, n) (left), and the associated eigenmodes 
(right) with n-dependent splittings, (b) Driven mode cor- 
relation function (7ia'(0) for thermal mechanical occupation 
A^th = 0, 1,2. Solid lines show the analytic calculation with 
7 — >■ 0, and dots show the full numerical results for TVth = 2 
only. The inset shows the minimal g'^aaiQ) as a function of 
A^th for several coupling strengths, (c) Same as (b) for the 
undriven mode correlation function (;ss'(0). Parameters are 
g/ hi — 20, and (for numerics) 7 = 0.001. 



IV. DELAYED COINCIDENCE AND SINGLE 
PHONON STATES 



In addition to the equal-time correlations discussed 
above, quantum signatures can also be manifested in pho- 
ton intensity correlations with a finite time delay. We 
now turn to a discussion of delayed coincidence charac- 
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terized by the two-time intensity correlations functions, 
"l(0)cl(r)c,(r)c.(0)) 



(36) 



for both driven and undriven modes, i ~ a, s. Express- 
ing this correlation in terms of a classical light intensity 
/, gC^){T) = (/(r)/(0))/(/)2, and using the Schwarz in- 
equality, we obtain the inequalities [46l |47] , 



\9^'Hr)-l\<W'H0) 



II 



(37) 
(38) 



Similar to the classical inequality 5^^' (0) > 1 at zero de- 
lay, violation of either of these inequalites at finite delay 
is a signature of quantum light. We calculate the delayed 
coincidence correlation functions for both the driven and 
undriven modes. 



A. Driven mode 

The correlation function g^a (t) is shown in Fig. [sj^a) 

for two values of the detuning A. The most striking fea- 

(2) 

ture is the apparent vanishing of gaa (t) at several values 
of T when the detuning is A = (curve A in Fig. pTa)). 
These are due to Rabi oscillations at frequency g/2 fol- 
lowing the detection of a photon. This vanishing of the 
finite delay correlation function is reminiscent of wave- 
function collapse that occurs in a cavity containing an 
atomic ensemble |46j . and while its origins are similar, 
there are important differences as we now discuss. 

We can understand the finite delay intensity correla- 
tions in terms of the simple six-level model discussed in 
the previous section. We extend this model to describe 
finite delay correlations by considering the effect of pho- 
todetection on the steady state of the system. Detection 
of a photon in the driven mode projects the system onto 
the conditional state ESI, 



|Ca|^)||' 



(39) 



where is given by Eq. (12) with steady state ampli- 
tudes and 1 1 • 1 1 denotes normalization after the jump. The 
conditional state \'ip°') has an increased amplitude Aioo 
after the jump (see jump at t = in Fig. [5] (b)). Follow- 
ing this initial photodetcction, the amplitude Aioo subse- 
quently undergoes Rabi oscillations with frequency g/2, 
and decays back to its steady state at rate 2k. For suffi- 
ciently large bunching at zero delay and strong coupling 
g > K, the Rabi oscillations of the amplitude Aiqq{t) can 
cause it to cross zero several times before it decays back 
to steady state. As the probability to detect a second 
photon is dominated by Aioo, its zeros are responsible 
for the zeros in the correlation function gaa {t) zero at 
these delay times. 




(b) 
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FIG. 5. 



glSir) for detunings A/g 



(a) Finite time delay intensity correlation function 
(A) and ^ (E). Detuning for 
curve A is the same as marked in Fig. [2] while E shows a 
new effect not seen at equal times. Thin dotted line indicates 
the bound (see Eq. (37l) for curve E. (b) Evolution of ampli- 
tude Aioo (normalized by its steady state value) at detuning 
A, A/g = 0, after detecting a driven Ca photon at r = 0. 
Vertical dashed lines mark delay times (ti, r2) where this am- 
plitude vanishes resulting in the vanishing of gSi (t) in (a). 
Parameters are g/n = 8 and 7/fi: = 0.02. 



The zeros in g'-^^r) appear similar to those exhib- 
ited in a cavity strongly coupled to an atomic ensem- 
ble [MlllTl [54] or a single atom 55J. However, in stark 
contrast to the atomic case, the zeros in Fig. [sja) are 
the result of Rabi oscillations following the initial quan- 
tum jump. This is qualitatively different from the atomic 
case, where the change in sign of the relevant amplitude 
(the analogy of Aiqq) occurs immediately after the jump 
itself, and the amplitude is damped back to steady state 
at the atomic decay rate F, without Rabi oscillation. As 
a consequence, the vanishing correlation function in the 
atomic case occurs at a delay set by tq ^ 7~^lnC, re- 
quiring only strong cooperativity C — g'^/Kj > 1 to be 
visible. On the other hand, the zeros in Fig. [sja) occur 
at delay times set by r ^ requiring strictly strong 
coupling g > K. 

Before moving on to correlations of the undriven mode, 
we briefly discuss the correlations of the driven mode at 
the other value of detuning shown in Fig. [5](a). At de- 
tuning A — ^ (curve E), which shows bunching at zero 
(2) 

time, gaa (0) ^ 1, increases above its initial value at fi- 
nite delay. This is a violation of the classical inequality 
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in Eq. (37), and is an example of "delayed bunching," or in Eq. (12|, 



an increased probability to detect a second photon at a 
finite delay time. A similar effect was recently studied in 
a single mode OMS [41 . However, like the Rabi oscilla- 
tions, the increased correlation function decays back to 
its steady state value of 1 on the timescale of 



B. Heralded single phonon states 

We now turn to a discussion of the delayed coinci- 
dence correlations of the undriven mode Cs- We note 
that correlations of the driven and undriven modes can 
be measured separately provided sufficient frequency res- 
olution, smaller than the mechanical frequency. The cor- 
relation function gls ij) of the undriven mode is shown 
in Fig. [6] for several values of detuning. Similar to the 
driven mode, the correlation function of the undriven 
mode exhibits Rabi oscillations that decay on the short 
optical timescale For detuning A = and IS./ g = ^ 

(curves A and D in Fig. |6|, the correlation {t) is de- 
scribed by our previous six-level model of Eqs. (13-18). 
However, at detuning A — 

„(2 



^ (curve E), we see that 

'(r) has a long tail that decays on the much longer 
mechanical timescale I/7. This is due to the heralded 
preparation of a single phonon by detection of a photon 
in the undriven mode, as we now discuss. 



^ CO 



A/5 = 
A/5 = 76/4 
A/5 = l/\/2 





D 



20 



40 60 



80 



100 



I V') - ■ ■ • + ^001 m)+ Awi 1 101) + A012 1012) , (40) 

since these are the states populated by detection of a Cs 
photon from the original six states (see Fig. [?]). Using 
the same approach as before, we obtain the following 
equations for the amplitudes, 



iioi 

^012 



7 



An 



001 



kAi 



(41) 

(42) 
(43) 



where we used j <^ k, and kept the leading term in 
Eq. (41). We obtain gss (t) by solving these equations 



for initial conditions determined by the conditional state 
I'lp'') after a quantum jump. 



(44) 



which is a superposition of states |001), |101), |012) (see 
Appendix for details), but in the limit of weak driving 
consist mainly of |001). 

Detection of a photon in the undriven mode implies 
that the three- wave mixing interaction converted a pho- 
ton from the driven mode into the undriven mode by si- 
multaneously adding a phonon. The relevant three-level 
subspace after the jump (see Fig. [7| has a similar struc- 
ture as in the steady state, but the presence of an extra 
phonon modifies the splitting of the one-photon states 
|1'±) = (|101) ± |012))/V2 to ^ (instead of g/2 with- 
out a phonon). This changes the one-photon resonance 




101 012 



V. 



steady state 
(off-rcsonant) 



conditional state 
(resonant) 



FIG. 6. Finite time delay intensity correlation function 
g^X-r) for detunings A/g = (A),V^/4 (D) and ^ (E). 
Thin dotted lines indicate the classical bounds (see Eq. ( |37[ )). 
Labels A, D, E correspond to the same detunings marked in 
Fig. [2]and Fig. [s] Parameters are g/n = 8, and 7/K = 0.02. 

The increase in delayed coincidence can be understood 
by extending the above analytic six-level model to ac- 
count for the conditional state of the system after detec- 
tion of a photon in the undriven mode. To do this, we 
simply add three additional states to the six-level ansatz 



FIG. 7. Effect of detection of a Cs photon at detuning E 
(A/g = ^). In steady state (gray region on left), the drive 
is far off-resonant. However, after detection of a Cs photon 
the system jumps into the conditional subspace (pink region 
on right). Due to the presence of an extra phonon in this 
subspace, the drive is resonant and the probability to detect 
a second photon is much higher than in steady state. This 
increased probability persists as long as the extra phonon, 
which decays slowly at rate 7. 

condition for the drive to A = Therefore, at this 
value of the detuning, the process of exciting the system 
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and emitting a single Cg photon is off- resonant; while af- 
ter the detection of a first Cs photon the system is pre- 
pared in 1 001), bringing it into resonance with the drive. 
This enhances the probability for subsequent excitation 
and emission of a second Cs photon, increasing the cor- 
relation function at finite delay. The maximum delayed 
coincidence occurs after a delay of r ~ 1/k, when the 
photons have reached the metastable steady state in the 
conditional subspace with one extra phonon. Eventu- 
ally, the delayed coincidence returns to its true steady 
state value of one on the timescale t ^ I/7, which is 
the mechanical decay time of the state |001). Note that 
the probability to detect a photon from the driven mode 
also increases in the conditional state, so a similar effect 
is seen in the delayed cross-correlation function gas (t) , 
where the photon in the undriven mode is detected first. 



V. REACHING STRONG COUPLING 

The nonclassical correlations predicted in this paper 
require strong optomechanical coupling, g > k, as well 
as sideband resolution, cj^ ^ k, 7. While the com- 
bination of these conditions has not yet been demon- 
strated, several experimental efforts are currently di- 
rected at reaching this regime. By using micro- and 
nano-fabricated OMSs such as microtoroids or photonic 
crystal beams, high frequency mechanical systems with 
ujm ~ 50 MHz - 5 GHz can be combined with low loss 
optical modes, such that the condition ujm S> k 3> 7 is 
satisfied [9, 2lJ At the same time the mechanical 

system can already be prepared close to the quantum 
ground state by working at cryogenic temperatures. In 
micro-fabricated OMSs, single-photon couplings of about 
g/K « 0.001 have been demonstrated [III |M1 ESj The 
largest value to date of g/n w 0.007 has been reached in 
photonic crystal beam resonators |56j , where colocalized 
optical and vibrational resonances are highly confined to 
maximize coupling while the surrounding structure is en- 
gineered to minimize loss. Conversely, in cold atomic ex- 
periments the effective strong OM coupling regime has 
been reached [57, , while sideband resolution remains a 
challenge [5H] . 

There are several existing proposals for how to meet 
the challenge oi g/n > 1 in the photonic crystal beam 
setup. First, the single-photon optomechanical coupling 
can be increased by making use of nanoslots in the struc- 
ture [m HO] to further localize the electric field at the po- 
sition of the mechanical mode. This could improve 5 by a 
factor of 10 [Hj. Second, numerical studies suggest that k 
can be further decreased by fine tuning the size and posi- 
tion of the slots in the photonic crystals (61] [62] . Finally, 
new materials are currently being tested for an overall 
improvement of the OM properties of nano-fabriquated 
devices [63]. Thus by using these ultrahigh Q photonic 
crystals or similar designs, an increase of g/K by a factor 
of ~ 100 is realistic. Note that once the strong coupling 
condition has been achieved, the implementation of two 



or multimode OMSs with adjustable tunneling 2 J ~ ujm 
can be realized via evanescent field coupling, as has al- 
ready been demonstrated in the weak coupling regime 

[211 [S3 HE]. 



VI. CONCLUSIONS 



We have studied nonclassical intensity correlations in 
a driven, near-resonant optomechanical system with one 
mechanical and two optical modes. In the regime of 
strong coupling g > k, this system allows for nonlinear 
quantum optics through a resonant three-mode interac- 
tion in which the exchange of two photons is mediated by 
a phonon. We have identified several different processes 
that can lead to nonclassical antibunching and delayed 
bunching, and we have derived a simple analytic model 
that allows us to describe and interpret photon-photon 
correlations in this system both at zero and at finite tem- 
perature. Our findings will be important as experiments 
approach the regime of strong OM coupling, and for po- 
tential applications of OMSs for quantum information 
processing. In particular, the long-lived correlation found 
for the undriven mode raises the intriguing possibility to 
exploit such a setup as a quantum memory. The genera- 
tion of heralded single phonons on detection of a photon 
from the undriven mode may have implications for build- 
ing OM quantum repeaters and quantum communication 
devices. 
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Appendix: Analytic model 

In this appendix we provide the analytic solutions used 
to calculate one- and two-time correlation functions in 
steady state. First, one-time correlations are calculated 
from the steady state solutions of Eqs. (13 18). We set 



the time derivatives to zero and solve the equations iter- 
atively, order by order in the weak drive. This procedure 
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yields 



^000 ~ Ij 
^100 = — 

Aqii — — 
^200 = ^ 

^111 

^022 = 



1 



2x 



1 + 4x2 ' 

1 + 

7i(l + 4x2)(l + 6a;2)' 

(1 + 4x2)(l + 6x2)' 
4x2 

(l + 4x2)(l + 6x2)' 



(A.l) 
(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 



where a = Vl/k (|a|2 ^ 1), x = g j (4k) and k ^ k — zA. 
Using these amphtudes, we can express all equal time 
averages. The mean photon numbers are 



na = l^iooP, 



Aioo + i 



(A.7) 
(A.8) 

(A.9) 



and the photon-photon correlation functions are 



2\A2 



noo\ 



i2)(n\ - ^l^022P 



1^011 1 



200 



+^ioo| 



(A.IO) 
(A.11) 

-. (A.12) 



To leading order in n/g these yield Eqs. (19 24 1. 

At finite temperature we calculate steady state am- 
plitudes within each phonon subspace n similarly in the 



ansatz of Eq. (26). Using the notation i?K(<^) introduced 



in Section IIIB the steady state amplitudes within the 



subspace with n phonons in the optical groundstate are 



\AiOn\' 



n^R^{g/V8) 



(gV(2n + l)/8) (f \/^^) 
r2V(" + !)(" + 2) 
32i?, (5v/(2n + l)/8) i?« (fVr^) 



(A.13) 
(A.14) 
, (A.15) 



|i20„P 



|^02n+2|^ 



(A.16) 

Two-time correlation functions are calculated sim- 
ilarly, using the conditional state after a jump (see 
Eqs. (391 and (44|) as the initial condition. For example 



the unnormalized state after detection of a photon in the 
ca mode is Ca\^) = iioolOOO) -)-\/2i2oo|100) +iiii|011). 
We solve Eqs. (13 T5| for the amplitudes with this state 
as initial condition. The finite delay correlation of the 
driven mode is 



9, 



(2)(^)^ \aMt)Y 



\A 



(A.17) 



100 1 



in good agreement with the numerics. The correlation of 

f 2) 

the undriven mode gls {t) is calculated similarly. The 
unnormalized state after detection in the Cs mode is 
Cs\-<p) =^ii|001)-h^ii|101)-hV2Ao22|012). Using this 
as the initial condition we solve Eqs. (|4l}[43| for the am- 
plitudes in the conditional state. In the limit of 7 ^ k 
we obtain 



l^oi2(r)|^ 



(A.18) 
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